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Abstract 

Problems of strongly interacting electrons can be greatly simplified 
by reducing them to effective quantum spin models. The initial step is 
renormalization of the Hamiltonian into a lower energy subspace. The 
positive and negative U Hubbard models are explicitly transformed into 
the Heisenberg and -x-xz models respectively. Basic tools of quantum 
magnetism are introduced and used: spin coherent states path integral, 
spin wave theory, and continuum theory of rotators. The last lecture con- 
cerns pseudospin approaches to superconductivity and superfluidity. The 
SO(3) rotator theory for the -x-xz model describes the charge density 
wave to superconductor transition for e.g. doped bismuthates. Analo- 
gously, Zhang's theory for collective modes of high T c cuprates describes 
the antiferromagnet to d-wave superconductor transition using SO (5) ro- 
tators. Finally, the Magnus force on two dimensional vortices and their 
momentum, are derived from the Berry phase of the spin path integral. 
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Part I 

Deriving the Effective 
Hamiltonian 

Let us consider the Hubbard model for conduction electrons hopping on a lattice 
with short range interactions 

n = r +u 

U = U^n, A n tA , (1) 

i 

It is always tempting to reduce the interaction term U to fermion bilinears 
(single electron terms) using the Hartree Fock (HF) variational approximation. 
However this approach is known to be seriously flawed in several important 
cases. For example, while the HF spin density wave is energetically favorable 
for U > at half filling, it breaks spin symmetry too readily in one and two 
dimensions, in violation of the Mcrmin Wagner theorem. This implies that Fock 
states might be too restrictive as a variational basis. We can illustrate this point 
using a simple toy model: the Hubbard model on two sites. It will also teach 
us something about onsite interactions and their effect on spin correlations. 



1 Two-site Hubbard model 

The two-site Hubbard model with two electrons, is 

H = -t^2 ( C M c 2,s + 4,s c m) + U n iA n hl> ( 2 ) 
s=T4 i=1 > 2 

where we take U > . The total spin S = Si + S2, commutes with the 

Hamiltonian, we restrict ourselves to the 5* = and 5* = 1 subspaces. The 
singlet subspace is spanned by three states 

|Q) ^ c ^ |Q) c y f ^ |Q) (3) 

In this subspace, the Hamiltonian is 



V2t Vzt 
H s =o = I V2t U 
V2t U 
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Figure 1: Noninteracting (U—Q) eigenstates of the two-site Hubbard model. 
The ground state is the "two electron Fermi sea" . 



Diagonalizing the Hamiltonian in the singlet sector, one gets 

/ (u - vu 2 + mt 2 ) /2 o o 

H s =o — ^ U _0 

V o o (u + vu 2 + mt 2 ) /2 

which, in the strong coupling limit (U 3> t) is 



H s=0 



u~s>t 



-tt 2 /U 
QUO 
U + 4t 2 /U 



The triplet subspace is spanned by the states: 



c i,T c 2a + c u c 2d| 0) ; 4 4 |o> , cl 4 4 a |o) 



V2 



(4) 



which all have zero energy. 

The comparison of the spectrum and wave functions of the non-interacting 
case and the strong coupling limit is displayed in Figs. |l| and S. 

The non-interacting groundstate is the two-site version of the Hartree-Fock 
two electron "Fermi sea". This state, and the lowest excited multiplet, contain 
sizeable contributions from the doubly occupied singlets 



v/2 



(IT4i,o 2 >, |0 1 ,T4 2 ». 



(5) 
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Figure 2: Eigenstates of the two-sites Hubbard model in the strong coupling 
limit. 



The Hubbard interaction pushes these states to energies of order U. At strong 
coupling therefore, the ground state becomes a valence bond singlet, of singly 
occupied sites with no charge fluctuations. It also cannot be expressed, even 
approximately, as a Fock state. 

The lessons to be learned from this toy model is that repulsive interactions 

can 

• enhance magnetic correlations, 

• reduce double occupancies in the ground state, 

• separate spin and charge excitations. 



2 Renormalization to Low Energy Subspace 

Low temperature phases and their interesting DC transport properties, are de- 
termined by the ground state and low energy excitations. This chapter is slightly 
formal, as it shows what is precisely meant by the Renormalization Group (RG) 
transformation which replaces an orginal non diagonal Hamiltonian with an ef- 
fective one for a lower energy subspace. Let us consider any Hamiltonian written 
as 



H = H + V , 



(6) 
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where Hq is diagonal and V is a non diagonal perturbation. We define the 
Hilbert space using the eigenstates of H , and Po(A) is projector onto the sub- 
space with energies less than A. 

The resolvent operator, G — (E — projected onto the latter subspace 

is given by a well known matrix inversion identity ffl 



G 00 (£) = P G(E)P = 

E-P (H + V)Pq 



Po 



Po 



PoV(l - P ) 

(l-P )VP Q (1-P )(E-(H + V))(1-P ) 
= {E-P (H q + V)P - 

- P V(1 - P )[(l - P )(E - (H + V))(l - Po)] ~\l - PoWPoY 1 

^(E-HejfiE))- 1 , (7) 

where the last equality defines the effective Hamiltonian Ti. e ff(E). Then the 
spectrum of H, which corresponds to states with non-zero weights in the sub- 
space considered, is given by the zeros of the characteristic polynomial of H. e f f{E), 



det (E n -n eff (E n )) 



(8) 



that is by the poles of the function TrG(£?)oo- The effective Hamiltonian can 
be also written as 



Heff(E) = P ^(Ho + V) + V(l-P ) 



(l-P )(l-(E-Ho)- 1 V(l-Po) 



(l-PoKE-Hoy^Po 



Po\(H + V) + Vj2 



1-Pa 
E-Ho 



V 



Po ■ (9) 



If Po projects onto the ground state manifold of Ho, Eq. ([|) defines the 
Brillouin- Wigner perturbation theory^. For two cut-off energies A' < A, Eq. (||) 
is a Renormalization Group transformation 



H{A) H(A') = Ho(A') + V(A') 



where 



H(A') = P(A')\ H(A) - V(A) J2 



n=l 



1 - P(AQ 
E-Ho(A) 



V(A) 



P(A') 



(10) 



(11) 



1 Note that the sum X!nLl m Eq. (|9|) does not correspond to a perturbation series for the 
ground state energy, since the terms depend on E. 
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After doing the best job we can to evaluate Ti(A') (it is clear one needs to 
truncate the infinite sum and do something about the energy dependence of the 
denominators), the terms separate naturally into a diagonal operator Hq and 
residual interactions V(A'). Sometimes we are lucky, and complicated terms of 
V(A') become relatively smaller as A' is reduced. These are called irrelevant 
interactions which scale to zero. We end this section by remarking that the RG 
transformation should preserve all the symmetries of the Hamitonian. If Ti has 
explicit symmetry-breaking terms, those terms may grow or shrink under the RG 
transformation rendering the low energy correlations less or more symmetrical, 
as the case may be. 



3 From Hubbard to t— J and Heisenberg Models 

As an explicit derivation of an effective Hamiltonian outlined in Sec. g, We 
consider the Hubbard model H = T+U of Eq. (||) in the strong coupling regime 
(U/t ^> 1). The diagonal part Hq we choose as U, the onsite interactions. This 
term divides the Fock space into two subspaces, the singly occupied and empty 
sites configurations 

S= I K iT ,n u ,n 2 , T ,n 2 ,x, . . •> /Vi,n iiT + n id < lj , (12) 

and configurations with one or more doubly occupied sites 

D = | \mj,n-L,i,n 2 ^,n 2 ,i, ■ ..) /3i,n iA +n iti = 2} . (13) 

The hopping term T couples the S and D subspaces by moving an electron into, 
or out of, a doubly occupied state. We define Pq to project onto the ground 
state manifold of subspace S, and thus 

Goo(E) - P Q G(E)P = (E — HeffiE))- 1 , (14) 

where the effective Hamiltonian H e f / , is given by Eq. (Q) 

n e ff(E) = p TP + p r[(i-p )(E-(u + r))(i-p )\ Vp . (15) 

In the strong coupling limit, expanding the effective Hamiltonian to zeroth order 
in E/U and to second order in t/U one gets 

H e f f (E) ► H l - J , (16) 

t/a<i 

n' 1 = p A T -jj E i^^T^ii-^ fo, (17) 

V i,j,k,s,s' I 
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i.e., the low energy excitations of the Hubbard model are described by the 
Hamiltonian of the so called t — J model. 

The fermion operators appearing in Eq. (|l7|), can be rearranged in the fol- 
lowing way: 

H*- J = P (T + T +H H ) P , (18) 



2 r^k 

r = --Y 

2U ^ 



c i,, c k,» n j - c\ac k ■ c\aCj 



(19) 



^(S.^-Sffl), ,20) 



(id) 

where J = 4t 2 /U and the S = 1/2 spin operators Si are 

V" i > • ( 21 ) 

a a being the Pauli matrices. 

At half filling, i.e., when rii = 1, Pq annihilates T and T since there can 
be no hopping processes within subspace S when there are no empty sites. The 
transport of charge is prevented by an energy gap of order U. This is the 
Mott insulator, which describes the undoped (parent compounds) of the high 
T e superconductors of the cuprate family. 

In this limit, the t-J model simply reduces to the spin 5=1/2 Heisenberg 
model 

H 1 - 1 -»■ - Si ■ Sj + const. (22) 

(id) 

As in the two site Hubbard model of two electrons (see Sec.|l|), the low energy 
excitations are purely magnetic. 



4 The Negative-U Hubbard Model 

The negative-U Hubbard model describes local attractive interactions between 
electrons which could be produced by several microscopic mechanisms e.g., 
phonons, plasmons or spin fluctuations. We choose, for convenience to write 
the model as follows 

w- u = -*£4 >sCj> -|$>-i) 2 + 

(id),s 1 

+ 1^2 V V ( Ui - - 1) - ni , (23) 

(id) « 
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where the negative-U term favors pairs of electrons on the same site in com- 
petition with the hopping term which delocalizes the electrons; Vij is intersite 
Coulomb interactions and /z the chemical potential. The following canonical 
transformation 



Ct,t — ► (k,-\ 
maps the negative-U to a positive-U Hamiltonian: 



(24) 

c i,i — 4a > (25) 



H -u ^ n+ u = _ t ^(ct T g. T _ + U - ^(n, - I) 2 + 

+ \ E J ii § ' § j -hY,S?-N(»+j), (26) 



where 



Jg = 4%, (27) 
h = 2fi, (28) 

= (29) 



and N is the total number of sites. 

Following the derivations of Sec. [| and using the fact that 7i +u is at half 
filling (for a proof see Sec. 3.3.1 in this model at large \U\/t can be directly 
mapped onto an effective pseudospin model 

h +u n~ x - xz + o(t 2 /u) 

1 nn 

W— " = -^[j^S|- J$(S?§f + §?§])] -£>S? ( 3 °) 

ij i 

where the pseudospin operators are 

St = \{ni-l) 



2 

Sf = ^{cl^+cucn) 



= ^"( 5 iT 5 a ~ 5 ^t) • ( 31 ) 

We see that the local charge operator and the pair operator have the same 
commutation relations as angular momenta along the z-axis and xy plane re- 
spectively. The quantum properties of the pseudospins explains Josephson com- 
mutation relation between charge and superconducting phase, [N, cf>] =1. 
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At weak coupling \U\/t < 1, it can be also argued that the negative-?/ 
model renormalizes onto a similar effective model as (|3^) albeit with different 
lattice constant and interaction parameters. The Fermi sea is unstable with 
respect to attractive interactions as seen diagrammatically by the divergence 
of the vertex function in the BCS or charge density wave channels. Since the 
attractive interaction scales to strong coupling, at the scale where the cut-off 
energy equals the BCS gap, the effective Hamiltonian can be transformed to 
the -x-xz model to obtain the strong coupling fixed point Hamiltonian. This 
procedure however, has not yet been carried out, to the best of our knowledge. 

In Lecture III, we shall use the classical ?{- x - xz model to describe super- 
conductivity, and charge density wave phases. 
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Part II 

Quantum Magnetism 

This lecture is technical in nature. It contains a brief review of the spin path 
integral and how to obtain the classical and semiclassical approximations to it. 
A fuller background for this subject, with compatible notations, can be found in 
Ref.Q. Here, a new emphasis is placed on anisotropic models and their rotator 
representation. 

5 Spin Coherent States 

Path integrals provide formal expressions which can lead to useful approxima- 
tion schemes. A path integral representation of spin models can be constructed 
using spin coherent states. Let us consider the eigenstates | S, m) of S 2 and 
S z with eigenvalues S(S + 1) and to, respectively. Spin coherent states are a 
family of spin states labelled by a unit vector f2 = (9,<fi), where 9 and <fi are 
the lattitude and longitude angles respectively. The are defined by applying the 
SU(2) rotation operator to the highest weight state^ in representation S: 

\n) s = e^ sz e iesv e-^ s "\S,S) 

,n^s v(S + m)\(S -to)! 

with 

u(0,(f>) = cos(6»/2)e i ^ 2 , (33) 
v(9,(f>) = sin (6>/2)e^ /2 . (34) 

Using Eq. (^2|), two useful identities can be readily proven: 
• The resolution of identity 

/l p27T +S 

dcos9 d(j)\n)(n\ = V \S,m) (S,m\ = I s ■ (35) 

1 J ° m=-S 



25 + 1 



4tt 



• The overlap of two states with closeby unit vectors 

(n[iV)~ + eC^-^w^')) (36) 

where 9 is the average lattitude of the two vectors. 



2 The phase e lSc ^ represents a gauge choice with one singularity on the sphere: at the 
south pole. 
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6 Spin Path Integral 

The partition function of a single spin with Hamiltonian H is 



Z = Tr [e- pH ] = Tr 



e-* H e-* H . 



<H 



N e =0/e times 



(37) 



with (3 — 1/T, T being the temperature. By inserting N e — 1 resolutions of the 
identity (|3^), and in the limit e — » 0, by expanding each exponential to first 
order, one gets 



N.-l 



Z~ dCl 1 ...dQ Nc Y[ (Q(T n )\l-eH\Cl(T n+1 )) 



(38) 



n=0 



with the boundary condition Oq = &N e - In the limit N e — > oo (and e — > 0, /3 
A^ e e = const.) Eq. (Bg) defines a path integral 



Z = It>Q(t) exp 



dr - cos 9(t))<1>(t) + H[Q(t)] 
The time dependent term in Eq. (j39|) 



(39) 



(40) 



derives from the overlap between coherent states (|3(j). It is known as the Berry 
phase of the spin history and it is geometric, i.e., depends on the trajectory of 
fl(r) on the unit sphere. In fact it measures the area enclosed by the path O(t) 
on the unit sphere (Fig. 0). 









iSw 


n 


= iS [ 






Jo 



The classical Hamiltonian is defined as 

H[(1(t)] = (Cl(r)\H\n(T)) 



(41) 



An implicit assumption in Eq. ( |39| ) is that the path integral is dominated by 
smooth (differentiable) paths. This turns out to be unjustified, since discon- 
tinues paths matter for the correct ordering of quantum operators. For that 
reason, path integral results should be checked whenever possible against oper- 
ator methods. Ordering ambiguities give rise to erroneous quantum corrections 
to energies and local spin correlations. They do not effect, however, long dis- 
tance and long timescale correlation functions. 

The partition function for a system of N spins is 



r N 

Z = & JJO^(r)exp 

i—l 



N 



dTH[{&i(T)}] 



(42) 
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For example, the nearest neighbor Heisenberg model partition function is 

r - i 

z H =j>vne^p zsVw[jV| — / Vfli-o. 

J ■ 1 2 Jo .. .. 



(43) 



The spin coherent states path integral (42) are convenient starting points for 
deriving semiclassical, i.e. large S, approximations. The integration variables 
are unit vectors, i.e. classical spins. The quantum effects enter through their 
time dependent fluctuations. Keeping JS 2 — > J fixed and sending S — > oo, 

suppresses the contributions of fluctuating paths with 57 ^ 0. This leaves in- 
tegration over frozen spin configurations precisely as in the classical partition 
function 

r N 

Z ► / TTdQie-' 3 *}] . (44) 

J i=l 

Now it is possible to use S as the control parameter for a systematic expansion 
of the partition function. In particular, applying the saddle point approximation 
(analytically continued to real time t — it) yields 



S 



6u> [f2j] 



dt 



SH 



= 



5Qi Jo SQi 
which are the Euler- Lagrange equations of motion for classical spins 

sdf (t) x t£(t) = 4-h \{iif} 



(45) 



(46) 
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7 Spin Wave Theory 



When spin symmetry is broken, either spontanously or by explicit symmetry 
breaking terms, it is quite natural to use the semiclassical expansion of the path 
integral. We consider small fluctuations around a classical spin configuration, 
Clf which minimizes if [f2]. 



To leading gaussian order, the partition function is approximated by 



Z ~ e-? H 
where 



"V 



VSn(r) exp 



N rP 
SV5 2 w[fij] - dT5 2 H[{£l t { T )}] 
„-_i Jo 



i=l 



N 



(47) 



(48) 



6 2 Hpi 



1 " fP „ X 

- / dr&i- Sfli x 5&i , 

2 i= i J o 

1 ^ - 5 2 H unA] ~ 



(id) 



(49) 



5Qi, which are perpendicular to 0, can be projected onto the two tangential 
unit vectors which defines the harmonic oscillator degrees of freedom 



qi = 5Qi ■ 4>i 
Pi = S5Qi ■ Oi , 

These variables can be used to represent the gaussian fluctuations as 

r/9 



(50) 
(51) 



-W\ J V Pi V Qi exp jT dr (iE^tZZ^ - 1 (q , p) ( * 



where if' 2 ' is a dynamical matrix of coupled harmonic oscillators 



where 



H (2) = 



Ki. 



M" 1 

%3 



K P 
P t M" 1 



d 2 H 



dqidq-j 
d 2 H 



q=p=0 



q=p=0 



(52) 
(53) 

(54) 
(55) 
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are the force constant and reciprocal mass matrices respectively and 

d 2 H 



dpidq. 



(56) 



q=p=0 



couples coordinates and momenta. Eq. ( J52| ) is the harmonic spin wave parti- 
tion function of any quantum spin Hamiltonian, whose classical ground state is 
known. By diagonalizing its action one readily obtains the spin wave excitation 
energies and wavef unctions, and spin correlations can be evaluated to the sub- 
leading order in The complexity of the calculation depends on the lattice 
symmetry of the classical ground state, i.e. such as the size of its magnetic unit 
cell. For example, for the Neel state, it is two lattice unit cells. 

For completeness, we work out the spin wave dispersion of the antiferromag- 
netic Heisenberg model with a Neel state given by 



_j (f,0) 1£ A 
(f.Tr) !£ B 



(57) 



where A and B are the two sublattices in which the lattice can be divided. The 
harmonic degrees of freedom are 



9i(t) 



<t>i{t) i e A 
7T + &(t) ieB 



and 



Pi(t) = S cos 9i{t) 
The dynamical matrix of the model fl43|) is 



(58) 



(59) 



j 



E 



ipi-Pjf 



s 2 



k, 9kj 



(<li -Qj) 2 = 






zJS 2 (l - 7k) 



9k 



where 



7k 



l^k-d 



(60) 



(61) 



z and d being respectively the coordination number and the vector connecting 
one site to its nearest-neighbours. The dispersion relation of the small fluctu- 
ations around the ground state configuration, i.e., of the spin waves, can be 
found solving the characteristic equation 

S' _2 zJ(l + 7k) 

zJ(l-Tk) 



det 



= 0, 



(62) 



and is Wk = ^zJy/l ~ 7k, for two distinct spin wave modes. 
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8 Continuum Theory for Anisotropic Models 

Spin wave theory is restricted to the ordered phases of the Heisenberg model. 
However, one can still use a semiclassical approach even in the absence of spon- 
taneously broken symmetry. A short range classical Hamiltonian is mostly sen- 
sitive to short-range correlations. Thus the Heisenberg Hamiltonian, in the large 
S limit, has at least short range antiferromagnetic order. In the path integral 
approach it is possible to utilize the short lengthscale correlations to define a 
continuum theory without assuming broken symmetry. 



In this section we spend some time preparing the ground for Lecture III. 
To that end, we derive the continuum theory for the anisotropic xxz model in 
a magnetic field. The resulting path integral will be later used in the context 
of quantum properties of superconductors. The continuum theory is shown to 
be equivalent to SO(3) quantum rotators. The rotator formulation is readily 
generalizable to SO(5) symmetry, which is the topic of section [l4|. Subsequently, 
for the isotropic case we review Haldane's mapping of the quantum Heisenberg 
antiferromagnet in d dimensions into the nonlinear sigma model (NLSM) in d+ 1 
dimensions, and the main results which can be obtained by that mapping. 

The first step is to parametrize the spins using two continuous vector fields 
h and L, 



= r)ih a (xi) 



\ 



L(xi 



S 



^ , (63) 



where r\i — e ln ' Xi has opposite signs on the two sublattices. Each pair of neigh- 
boring spins (4 degrees of freedom) is replaced by n, L. We can choose n to be 
a unimodular (|n| 2 = 1) Neel field (2 degrees of freedom), and L is the perpen- 
dicular canting field, with the constraint L ■ h = (2 degrees of freedom). The 
spin measure of Eq. (|42|) becomes 



VVL -> VhVLS 



(£ ■ nj (64) 



where the 5 functional are local space-time constraints. 

Let us consider a general anisotropic spin model in a magnetic field h, 



n = \Y J JffiCi«-h-sY J ^. (65) 



2 



Using ( |63|) and expanding to lowest order in L a , din a and diL a we obtain the 
energy density 

TL — > J d d x H \Z, nj 

H = E cl [n] + ( L ° f + P^ n °) 2 ) a ~ d ^L, (66) 
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where the energy, spin stiffness, and susceptibility parameters are respectively: 

= V j° 



'tjWl^i-^jl 2 ) 



a and A/ - are the lattice constant and size respectively. Expansion of the Berry 
phase term to the same order yields two terms 

iS^w^ = -iT + ij dr J d d x nxh- L, (68) 

where 

T(n) = 5 /" (fW^^z)) . (69) 
Collecting the terms together, we have the path integral 

Z xxz = J VnVLS (L ■ nj e~ lT cxp (J drd d x in x h ■ L - H[L, n]j . (70) 

9 Anisotropic Quantum Rotators 

Eq. (|7^) can be physically understood as a path integral of rotators. Consider 
the phase space path integral over an N dimensional field ft, and canonical 
momenta p with a "Mexican hat" potential 

, - J WDR exp (J drd^x iit-ff- H\p,n]j 

h mh = H rot \p,n]+K{\n\-lf . (71) 

If K is taken to be very large, fluctuations of drin = \n\ — 1, and its conjugate 
momentum p\i become high frequency harmonic oscillators, which can be inte- 
grated out in the adiabatic approximation. This leaves us with the slow degrees 
of freedom n — n/\n\, and p±, and a renormalized Hamiltonian 

Zk^oo oc J VhVpS(p- n) exp (J drd d x in ■ p ± - H rot [p ± ,n^j . (72) 
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A Faddeev-Jackiw quantization of a particle on an ./V-sphere leads to the same 
constraints, as shown in Ref. ||. 

For the N=3 model, the angular momenta and the transverse momenta are 
both vectors, related by 

L = n x px, p± = L x n. (73) 

A change of variables p± — ► L has a unit Jacobian 

VpS(p- n) = VLS{L ■ n) . (74) 

Substituting @ and (f74|) into |7|) yields a path integral of the form @, 
without the Berry phase e~ lT . 

For general N, the N(N — l)/2 angular momenta are defined as 

L ab = n a p b -n b p ai a,b=l,...,N. (75) 

Note, that it is not useful for N > 3 to write the path integral measure in 
terms of L abl because they are not independent degrees of freedom, and more 
constraints are required. On the other hand, the transverse momenta, which 
obey Yl a n a(p±)a = 0, can be expressed as 

(p±) a = ^L ab n b . (76) 

b 

Therefore it is always possible to express the rotator Hamiltonian H rot of ([72]) 
in terms of L ab and h. 

The opposite direction, might be even more useful. For example, as in Eq. 
( p6[ ) , the starting point could be a Hamiltonian whose kinetic energy is expressed 
using symmetry generators 



H r °%h] = E*[n] + lj2x- b 1 L 2 ab + \j2p^ n ^ 

a<b a 

- d Y, h abL ab , (77) 



a 

a<b 



where h ab and \ab are SO(N) fields and susceptibilities respectively. (For SO (3) 
their vector notation is given by X a = ~}2 bc ^ abc X bc ). Substituting ( |75| ) into 
H rot yields 



H rot [p,h] - E cl [n] +±J2 M abPaPb 

ab 

+ Ps (^ n ") 2 _ a <1 X! h ab(n a Pb ~ n b p a ), 

a a<b 

M- b l = 6 ab (^2xacn 2 )j -X~ b l n a n b (78) 
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where M[n] is an anisotropic "mass " matrix in the Cartesian basis. 

We note that by (|T8|), the path integral (|7^) is Gaussian in momenta p a . 
One must be careful and integrate only over the transverse components to n. 
For a given direction n, we choose the transverse basis ej, i = 1,N — 1 which 
obeys the following conditions 

ei-n = 0, i = l,iV-l 

N 

Y,% M ob% = ^M^[n]. (79) 

This is always possible since the first condition leaves the freedom to perform an 
SO(N-l) rotation on the transverse basis. For an arbitrary transverse basis {/»}, 
we find the rotation which diagonalizes fiM~ 1 fj, and the resulting eigenbasis 
is chosen as {e^}. 

Thus, we parametrize 

f=J2PiCi (80) 

i 

and integrate unrestrictedly over J Vpi to obtain 

Z = I Vncxp (- f dr f d d x(E cl [n]+^-p a (din a f 
J \ Jo J ; a 2 

+ \Y. M i -^)) 2 )) 

3 / 

hj = a- d J2 h ab(n b e] - n a e)) . (81) 

a<b 

This expression is ready for the evaluation of the classical (time independent) 
ground state fi cl as a function of applied field h a b, and a spinwave expansion 
about it. In the next section however, we shall only deal with the isotropic case. 



10 Haldane's Mapping 

We now return to ( |Sl| ) but specialize to the isotropic Heisenberg model without 
a field, where life simplifies considerably. 
The inverse mass matrix is simply 

M~ l = x'^ab - n a n b . (82) 

Any choice of transverse basis ii yields a diagonal Mj = x in Eq. (ff^) . Omitting 
the constant E cl , and inserting the Berry phase e~ lT ^ the partition function 
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( |Sl| ) reduces to Haldane's result 

Z = J Vhe-^cxp (y\[ dT J d d x( X \M 2 + PsJ2\ d ^\ 2 ^ = 

= J Vhe- ir ^exp f-^r- J dx d d x {d^hf^j , (83) 



where xq = ct, and c = yps/x is the spin wave velocity. For the nearest 
neighbor model with interaction J, one obtains from (|6^): 



X = S 2 /(MJa d ) , (84) 
PS = Ja 2 - d , (85) 

/ = ^— = 2V2S- 1 . (86) 

PS 

Eq. (|8^) is the partition function for a NLSM with an additional Berry phase 
term. 

In 1 + 1 dimensions the NLSM is disordered for all /, as required by the 
classical Mermin Wagner theorem. Its correlations are known to fall off expo- 
nentially at large distances with a correlation length which goes as £ oc e 27r ^. 
By the (Lorentz) symmetry of the action between spatial and temporal dimen- 
sions, this implies a gap (Haldane's gap) for all excitations above the ground 
state. However one should also consider the effects of the phases brought about 
by the term T(n). 

For d = 1, T(n) is a topological winding number of the two dimensional 
NLSM. For all continuous fields, it yields e~ lT ^ = e~ t27:Sk with k an integer 
number. Thus, the Berry phase factor is 1 for all integer S, while it can be 
±1 for half integer spins. As a result, it produces interference effects for half 
odd integer spins, and drastically changes the ground state properties and the 
elementary excitations spectrum of the Heisenberg chain. 

In d = 2 the topological phase is zero for all continuous fields. For the nearest 
neighbours Heisenberg model, Neves and Perez || proved that the ground state 
is ordered for all S > 1. Also, series expansions and numerical simulations 
provide evidence of the presence of an ordered ground state even for S = 1/2. 



11 Spin Liquid States 

For most antifcrromagnetic Heisenberg models the ground state is not explicitly 
known. While for finite bipartite lattices, ground state theorems require the 
ground state to be a total singlet, and have positivity conditions (Marshall's 
signs). The discussion of the previous section expects them to exhibit long 
range order in two dimensions. Considering these conditions, particularly useful 
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variational states for the S = i antiferromagnetic Heisenberg model are Fazekas 
and Anderson's Resonating Valence Bond (RVB) states || 

\{d a }) = $> Q |a) , (87) 

a 

where 



ieAJeB 




and d a can be chosen to have the form 

d a = Y[ u i:j . (89) 

\(i,j)eA a J 

Since the states \a) are not orthogonal to each other, it is not possible to 
evaluate correlations of RVB states analytically. Monte Carlo simulations of 
Liang, Doucot and Anderson ||, and HavilioQ, have found that the RVB states 
have long range Neel order for iiy that decay slower than 

Uij ~ \xi — Xj\~ p , p < 3 . (90) 

The RVB states Eq.(|87|) can thus be used as variational ground states for both 
ordered and disordered phases. This makes them appealing candidates for 
studying the transitions from the Neel phase to possible quantum disordered 
phases, particularly in the presence of hole doping^. 
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Part III 

Pseudospins and 
Superconductivity 

The -x-xz model, encountered in Lecture 0, will be considered as an effective 
Hamiltonian for the low temperature, long wavelength properties of s-wave su- 
perconductors and charge density waves. Such transitions are observed in e.g. 
doped bismuthates Bai-^K^BiBOs. Recall Eq. (f{0|): 

nn 

n-'~" = jE[w-J§W + ^)]-E^- (si) 

ij i 

S are spin 1/2 operators, and h — 2/i, where /i is the electron chemical potential. 
As shown later, it can also be used to describe superfluids. In this lecture we 
use the continuum rotator model, and its classical limit to obtain the phase 
diagram. The generalization to Zhang's SO(5) rotators, which describe the 
collective modes of high T c cuprates, and the transition from antiferromagnetism 
to d-wave superconductivity, is straightforward. The last section uses the -x-xz 
model to describe dynamics of vortices in superfluids. 



12 Charge Density Wave to Superconductivity 

In bipartite lattice with sublattices A and B (e.g. square, simple cubic, etc.), the 
negative J x terms of the -x-xz model can be rotated to positive terms by a global 
e iS z TT ro t a ti on on sublattice B. Obviously, for frustrated (non bipartite) lattices 
the xxz and the -x-xz models are not equivalent. For example, the triangular, 
Kagome, and face centered cubic lattices Q , the -x-xz model prefers to order in 
the xy plane even when J 2 > \J X \ at zero doping^. 

Here we specialize to the classical model on a bipartite (square or cubic) 
lattice 



H = 



Hn 



Z (i,3) 

-hsJ2&! +H nnn 

i 

1 nnn 



z' 

* ,3 



(92) 

3 This helps to explain superconductivity in K3C6O which is an FCC compound with a half 
filled conduction band. 
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Figure 4: Typical mean field phase diagram (doping concentration versus tem- 
perature) of the -x-xz model. The arrows denote directions of rotator field fi in 
the xz plane. 

where z and z are nearest neighbor (rm) and next nearest neighbor (nnn) 
coordination numbers respectively. 

By ([H]) we note that the electron charge expectation value and the super- 
conducting order parameters are 

(ni)-l = n? 

(4Ai> = nf + in? 



x = -N-^&l (93) 

x is the hole doping concentration away from half filling. 

At finite temperatures, molecular mean field theory for S = 1/2 |J gives the 
critical temperatures as a function of doping x to be 

T? DW {x) = \(J*-K)(l-x*) 

T? c (x) = - JXX ^ (94) 
c w 4 arctanh(x) v ; 

The two curves meet at T* (See Fig. ||[). This point is tetracritical or bicritical, 
depending on whether the transition (as a function of field h) is second or first 
order respectivley. In Fig. |] this is reflected by the nature of the intermediate 
region which would be either a mixed ( "supersolid" ) phase, or phase separation 
between pure CDW and SC domains respectively. The next section is devoted 
to determining the criteria for the order of the transition. 

Near T* , the lowest order expansion of the Ginzburg Landau free energy 
functional has SO (3) symmetry. Thus even though the model might have high 
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anisotropy J z >> J x , near the bicritical point the Heisenberg symmetry is 
approximately restored. This "symmetry restoring" is also argued to happen 
near the multicritical point of the anisotropic SO(5) model pp[|. 

13 Order of Transition from Rotator Theory 

Finding out the order of the CDW-SC transition at zero T, requires the energy 
of the putative mixed state (M), which interpolates between the pure CDW and 
pure SC. The rotator partition function ( jgj| ) comes in handy for that purpose. 
We write ( |92| ) as ( |65| ) by letting — J x — > +J X , and take the continuum limit 
with unit lattice constant a = 1 following the perscriptions of section ||. The 
M ground state is parametrized by one Neel angle n = (sin 8, cos 8), where 
< 8 < tt/2. 8 = 0, = vr/2 are the pure CDW, SC respectively. 
First we evaluate the inverse mass matrix of (fTq ) 

Xy 1 cos 2 8 — Xy 1 cos # sin f 

M~ b 1 =[ x^cos^ + x^sin 2 ^ " 

-Xy 1 cos sin Xy 1 sm2 9 

(95) 



The transverse vectors which diagonalize the projected inverse mass matrix are 
ei = yx h and e.2 = y- The mass eigenvalues are 

" I= { X (x-cos^ + x-sin^)-)' (96) 



Set the lattice constant to unity a = 1, using (67) and (B2J) we obtain 



hj = —h sin 8 8j t 2 
X~\0) - ^((J z + K) + {J Z -K) cos 2 0+\J x \sm 2 0) 

Xy\0) = ^(\J x \ + (J z -K)cos 2 8+\J x \sin 2 8) . (97) 

Here we have set the lattice constant to unity a = 1. 

Thus, the classical ground state 8 cl minimizes the energy 



E rot [8, h] = - (J z - K) cos 2 9-J x sin 2 8 - 1 ^ M^ 2 



= -(J"-X) + (J Z -X- J x )shr 



3 = 1-2 

\S 2 h 2 sm 2 8 
2 a 4 



J z - K + J x + 2K sin 



2 , 



(98) 

A simple analysis of ([)8|) reveals that the order of the transition depends on the 
sign of K: 
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For K < 0, 

d 2 E 



<9(sin 2 <9) 2 



< 0, < 9 < tt/2 (99) 

K<0 



i.e., there is a first order transition between cl — to 6 cl — tt/2 called a "spin 
flop". This happens at a magnetic field /i sc given by 

E ro %h sc ] = ^[vr/2,^] 



-> /i sc = 25" V( JX - K - J X )(J Z + K + J X ). (100) 
For any h, the doping concentration is given by 

*(*) = 25- (101) 

In the case of a first order transition, there is phase separation between the 
doping concentrations of pure CDW and pure SC at densities 



/ J z — K — J x 

x cdw = 0, x sc = y jz - - - jx . (102) 

On the other hand, for K > 0, E(9) is minimized at 



sin 2 6 d = i^hS^I j z _ K K + _ J j x ~(J Z ~K + ,nj /2K (103) 

which indicates the existence of a mixed phase for fields in the range h c( i w < 
h < h sc , where 

2 

h cdw = gV(J z - K + J X ){J Z - K - J*j, x cdw = 0. 



2 / J z — K — J x J z — K — J x 

k - = S tp-K + J x ( JZ + K + JX )> x - = \j jz-K + J* - 

(104) 

For the pure nearest neighbor model [K = 0), it is easy to see by (|9^), that 
at h sc = h c d w the mixed state is degenerate with a Maxwell construction of 
phase separation into SC and CDW. This degeneracy is lifted in the quantum 
version of the same model, which favors phase separation |ll). 

14 SO (5) Rotator Theory and High-T c Super- 
conductors 

Recently an SO(5) theory of high T c cuprate superconductors has been proposed 
by Shou-Cheng Zhang |HJ . The order parameters of antiferromagnetism (AFM) 
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and d-wave superconductivity (dSC) are written as the cartesian components of 
a 5 dimensional vector 



n = 
(n 2 ,n 3 ,n 4 ) = 

rix + in§ = 



(711,712,713,714,715), 



1 



(105) 



g(p) is the Fourier transform of the pair wavefunction which for d-wave pairing 
transforms as cos(p x )— cos(jp v ) under lattice rotations. Zhang has found explicit, 
second quantized constructions for the 10 SO(5) generators {L a ^} a b, a,b = 
1, . . . ,5 which rotate n in 5 dimensions. Particularly interesting are the II 
operators, i.e. generators Li, a ,a — 2,3,4, which rotate between the AFM and 
dSC hyperplanes. These are expected to create new low lying pseudo-Goldstone 
modes near the transition. 

Zhang has proposed that the long wave fluctuations are described by an 
effective rotator Hamiltonian of the form M, 



H[L,n] 



1 



2 ^ 

a<b 



X ab 



1 



py,\^ri a \ 



2 r ^ 

a 

,2 1 „2 



+ 5(711 +n 5 ) - 2juLi 5 



(106) 



where L a {, are the generators of SO (5) algebra, and L15 is the charge operator 
whose expectation value yields half the doping concentration (£15) = x/2. Us- 
ing the substitution (75), the momenta can be integrated out of the partition 
function, leaving us with 



Z 



T>h exp 



~Io dr I d Hl^ Mj ( h " 

ip^|Vn Q | 2 +.g(n? + n2))j 

a=l / 



— ihj 



2/i(n 5 e) - nief ) 



(107) 



We set h = 0, and search for the classical ground state. 

Let us first consider the SO(5) symmetric model, where all Xab — Xi an( i 
g = 0. For any finite /x ^ 0, the ground state flops into (711,715) plane, i.e. is 
superconducting. 

Experimentally, a transition from AFM to dSC is observed at low hole 
concentrations x in many of the high T c systems. Appealing to the anal- 
ogy with the -x-xz model, this suggests that the symmetry breaking terms of 
SO(5) — > SO(3) x 5*0(2) should not be forbiddingly large. Symmetry breaking 



15 VORTEX DYNAMICS IN SUPERFLUIDS 



27 



terms can be included by g > 0, and letting the charge susceptibility \ c = xi.5 
be different than all other susceptibilities Xab — X> a >b 1, 5. Without loss of 
generality, we can choose the order parameter to tilt between the AFM and dSC 
hyperplanes n = (sin 9, cos 0, 0, 0). Following the same derivation as for xxz 
rotators (|9S|), the classical energy is 

E[8, /i] = g sin 9 2 ^fj^ . 2fl ■ (108) 

X~ cos^ u + Xc sin 

It is now straightforward to verify that the ground state is in the (712,713,714) 
sphere at /i = 0, and will "flop" into the SC state 8 — tt/2 at large enough [i. 

It also follows, using the same path as in Section [l|, that the order of the 
transition depends on the relative magnitudes of susceptibilities^: 



Xc > X =^ First order transition at /i^sc — \J g/(^Xc) 

phase separation at < x < 
Xc<X Mixed phase for y/g/(2 X ) < n < Vffx/S (Xc 1 + X^ 1 ) 

at densities < x < 4 a/ gx/2 

(109) 

In the mixed phase, the relation between the SC order parameter and the doping 
concentration is linear 

(nl+nl) = sin 2 6 cl = —^= (110) 

Before attempting to compare these results to experiments, we must remember 
that this is merely the classical approximation. 

15 Vortex Dynamics in Superfluids 

We heneceforth restrict ourselves to two dimensions. The quantum -x-xz Hamil- 
tonian ( |30| ) for spin 1/2 can be written in terms of Holstein-Primakoff bosons 
defined as 

S z = Ja-±, (111) 



5" = y/l-ataa , (112) 

S + =aWl- a f a , (113) 
4 These results differ somewhat from the phase boundaries in Ref. |lo| 
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yielding 

- JZ[a\a jy /(l - a\ ai ){\ - a]a 3 ) + h.c] - - a\ ai ) . (114) 

i 

The partition function can be written as a Bose coherent state path integral 



Z = J X> 2 zexp -i J dr ^ z*5 r z, - H~ 



x — xz r„* 



(115) 



where 2?Zi is an integration over the complex plane. Taking the continuum limit, 

Zi — ► 4>(xi), we can write 



Z = J P 2 0e _s[01 



where 



S = / drd d x 



0^ T 0+-a|V0| 2 + y(H) 



(116) 



(117) 



is the time dependent Ginzburg-Landau action with 

y(<^*)=zjvQ-|0| 2 ) -zjV|0| 2 (l-|^| 2 )+ M |0| 2 . (118) 

The classical equation of motion is given by analytically continuing r — > it and 
finding the saddle point 



= -> = aV 2 + ^ ) , 



(119) 



which is known as the Gross-Pitaevski or Non Linear Schrodinger equation, 
whose solutions 4>{x,t) describe collective modes (phase fluctuations), and dy- 
namics of vortex configurations. 

A di ffere nt approach to dynamics of superfluids, is to use the quantum spin 
model (114) on a lattice with a lattice spacing a which is smaller than the 
interparticle spacing, as represented by the spin coherent states path integral 
©■ 



Z = / T> cos 9iT>(f>i exp 



f ' dT M^(l-cos^(r))^(r)-iJ[a] 



(120) 
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Figure 5: Berry phases due to a moving vortex. 



A superfluid state is characterized by an ordered state in the xy plane, with 
a constant two dimensional boson number density a 2 p s = |0| 2 . Using (|l!3| ), the 
average spin direction is related to p s by 

(cos0)) = Ps a 2 . (121) 

A vortex configuration can be parametrized by the azimuthal angles at lattice 
points i by 

&(i) = argfo-avW) (122) 

where xy(t) is the vortex core trajectory. As one can see in Fig.|^, the Berry 
phase of a vortex path xy (t) can be written 




(123) 



where the sum is extended to the N c lattice sites included by the vortex path 
since the contribution of the others is zero. This Berry phase generates a Magnus 
force on the moving vortex. This is evident when we write it as an integral over 
a gauge potential 

dxA -'x = A- dl = BS C (124) 
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where S c is the area included by the path c, and B is an effective magnetic field 



(in dimensions of unit flux quantum) . Thus, comparing Eq.(123) and Eq.(124) 
one can define B to be simply 



B = Ps 

and the Magnus force acting on the vortex can be written 

FMagnus — Bz X Xy — 2np s Z X Xy . 



(125) 



(126) 



In the absence of any other time derivative terms, the vortex moves like a 
massless particle in a strong magnetic field restricted to the lowest Landau 
level. The semiclassical momentum p of a vortex configuration can be evaluated 
by computing the expectation value of the translation operator Tg 



■ p _ 



(Cl\Tg\Cl) = (Cl{x)\Cl(x + a)) 



exp 



exp 



iS ^(1 - cos 9i) (4>(xi) - <p(xi + a)) 



d 2 xpsW4> ■ a 



gi27r p s a- zxxv 



(127) 



Similarly, the momentum of a vortex-antivortex pair at positions xy, xy can be 
computed 

P = 2np s z x (x v - x v ) . (128) 
Since on a lattice the only distinguishable momenta are within the first Brilluoin 



zone, Eq. (128) implies that vortex-antivortex pair configurations can tunnel 
between different separations which belong to the discrete family 



x"y=-(2irp s 



P + G) XZ, 



(129) 



where G is any reciprocal lattice vector. This is precisely an Umklapp scattering 
of the superfluid current by the lattice. This amounts to quantum dissipation 
of the supercurrent due to continuos translation symmetry breaking of a lattice 
potential Q. 
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